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A NEW SUBSPECIES OF POCKET GOPHER (GEOMYS) FROM 


TEXAS (MAMMALIA: RODENTIA: GEOMYIDAE) 
Michael J. Smolen, Richard M. Pitts, and John W. Bickham 


Abstract. —' Two isolated populations of pocket gophers were discovered far 
outside the previously known distribution ofgeomyids in southern Texas. These 
represent a new taxon, which is related to Geomys texensis in chromosome 
morphology and biochemical properties. The two share identical diploid num- 
ber (2N = 70), fundamental number (FN = 68), and the presence ofa distinctive 
large acrocentric X-chromosome. There are no fixed differences in any of the 
18 loci analyzed with starch gel electrophoresis, when comparing the new taxon 
with G. texensis. However, both of these taxa share alleles at two loci that are 
present as fixed differences when compared to G. bursarius major. Analysis of 
cranial morphology indicates that the new taxon resembles G. texensis, but 
differs significantly in cranial dimensions of size and shape, related primarily 
to measurements of basal and palatal lengths and mastoid breadth. Because of 
the morphological distinctiveness and the extreme spatial separation (120 km) 
from populations of G. texensis confined to the central basin of the Edwards 
Plateau, we conclude that this new taxon is an isolated, relictual population of 
G. texensis, and is a distinctive subspecies. 


During a recent survey of the pocket go- 
phers in southern Texas, populations were 
discovered far outside the previously known 
distributional ranges of Geomys. Geomys 
attwateri and G. personatus are both broad- 
ly distributed species (Fig. 1) that occur over 
much of the eastern and central portions of 
an area herein referred to as South Texas 
and bounded by the Edwards Plateau, San 
Antonio River, Gulf of Mexico, and Rio 
Grande River (Davis 1940; Kennerly 1954, 
1959; Williams & Genoways 1981). Their 
distributions are restricted to sandy and 
sandy-loam soils (Davis 1940, Honeycutt & 
Schmidly 1979) that are dispersed through- 
out this region. The northern and western 
parts of South Texas, however, consist al- 
most entirely of hard indurate soils, with 
high clay content, and caliche. Geomys is 
occasionally found in isolated pockets of 
sandy or sandy-loam soils in this area. 

Survey of these isolated patches of suit- 
able habitats in the northwestern part of the 


region detected previously unknown pop- 
ulations in Medina, Zavala and Uvalde 
counties. Karyotypic, morphometric and 
electrophoretic analyses identify these as re- 
lictual populations that are related to Geo- 
mys texensis, and they are described herein 
as representing a new subspecies. 


Materials and Methods 


Animals were live-trapped (Baker & Wil- 
liams 1972) and returned to the laboratory. 
Metaphase spreads were prepared from bone 
marrow cells, using 0.075 M potassium 
chloride hypotonic and Carnoy's fixative 
(Baker et al. 1982). Fresh slides were pre- 
pared by dropping (1.5 m) the cells onto 
slides flooded with distilled water. Standard 
karyotypes were analyzed following staining 
with 296 giemsa stain in 0.01 M phosphate 
buffer. Diploid number (2N) was deter- 
mined by counting at least 10 spreads, 
whereas fundamental number and chro- 


6 PROCEEDINGS OF THE BIOLOGICAL SOCIETY OF WASHINGTON 


102 101 100 99 98 97 


ы n — T 
= { 


7 
И 


4 jp 
n | | 
| 


p 
| | 


ТИ 
MINA 
Г 


| 


28 A G. personatus 
0 100 
hee 
km 
27 
102 101 100 99 93 


Fig. 1. Map of the distributions of pocket gophers in southern Texas. Open circles represent localities from 
which specimens were pooled to form samples. 
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mosome morphology were described from 
photographic prints of selected spreads. 

Samples representing the unknown taxon 
(n = 21), G. texensis (n = 13), and С. bur- 
sarius (n = 10) were used to assay biochem- 
1cal variation. Heart and kidney tissues were 
minced in a grinding solution (Tris/EDTA/ 
NADP) and homogenized using a mechan- 
ical homogenizer (Tissue Tearor, Biospec 
Products). Samples were loaded into 1296 
starch gels (Starch Art). The techniques used 
for visualizing the allozymes were those de- 
scribed by Harris & Hopkinson (1976) and 
Honeycutt & Williams (1982). Proteins were 
examined on the following buffer systems: 
Poulik, continuous Tris-citrate II (pH 7.0 
and 8.0), Tris-malate EDTA (pH 7.4), and 
Ridgway (pH 6.7). Eighteen presumptive 
loci and their respective Enzyme Commis- 
sion numbers (Murphy et al. 1990) are as 
follows: glycerol-3-phosphate dehydroge- 
nase (G3PDH; Е.С. 1.1.1.8), malate dehy- 
drogenase (MDH- 1,2; E.C. 1.1.1.37), 1soci- 
trate dehydrogenase (IDH-1,2; E.C. 
1.1.1.42), phosphogluconate dehydroge- 
nase (6-PGDH; Е.С. 1.1.1.44), superoxide 
dismutase (SOD-1,2; E.C. 1.15.1.1), purine- 
nucleoside phosphorylase (PNP; E.C. 
2.4.2.1), aspartate aminotransferase (AAT- 
1,2; Е.С. 2.6.1.1), creatine kinase (CK; Е.С. 
2.7.3.2), phosphoglucomutase (PGM- 1; Е.С. 
5.4.2.2), esterase (EST-1,2; Е.С. 3.1.1.-), 
peptidase (РЕР-1 leucyl glycyl glycine, 
PEP-2 leucyl alanine; E.C. 3.4.-.-), aconitase 
hydratase (ACOH; E.C. 4.2.1.3), glucose- 
phosphate isomerase (ОРТ, Е.С. 5.3.1.9). 
The allozyme data were scored in a side by 
side comparison of mobility where identical 
mobilities were recorded as individual al- 
leles. Data were analyzed using BIOSYS-1 
(Swofford & Selander 1981), which pro- 
duced Rogers' genetic similarity and genetic 
distance values for the three taxa (Rogers 
1972). 

The new taxon was compared to G. attwa- 
teri, G. personatus, and G. texensis using 
univariate and multivariate statistical anal- 


yses of morphological data. Geomys attwa- 
teri was represented by two populations from 
the western extent of its distribution and 
thus in relatively close proximity to the new 
taxon. Geomys personatus streckeri was se- 
lected as a representative of the personatus 
complex because its distribution also ap- 
proaches that of the new gopher. Two pop- 
ulations of G. texensis represent the previ- 
ously recognized subspecies, G. bursarius 
texensis and G. bursarius llanensis. They 
were recently elevated to specific ranking 
and synonymized as the monotypic G. tex- 
ensis (Block & Zimmerman 1991) based on 
patterns of biochemical variation. Although 
neither G. texensis nor G. bursarius are 
closely distributed to the new taxon, they 
are included because they are karyotypically 
identical to it. 

Geomyids show extensive variation in 
secondary sexual characteristics (Baker & 
Genoways 1975; Honeycutt & Schmidly 
1979; Williams & Genoways 1977, 1978, 
1980, 1981), which necessitated separate 
analysis of males and females. Three exter- 
nal measurements were recorded from spec- 
imen labels, and 11 cranial characters were 
measured using dial calipers. Skull mea- 
surements were taken as described by Wil- 
liams & Genoways (1977). Measurements 
we evaluated are as follows: total length (TL), 
length of tail (T), length of hind foot (HF), 
condylobasal length (CBL), basal length 
(BL), palatal length (PL), prefrontal depth 
(PFD), length of nasals (LN), diastema 
(DIAS), zygomatic breadth (ZB), mastoid 
breadth (MB), squamosal breadth (SB), ros- 
tral breadth (RB), interorbital constriction 
(IOC), breadth across third molars (M3), 
length of maxillary toothrow (M1). Only 
adult animals were used and they were iden- 
tified by the degree of ossification of the 
sutures of the skull, particularly the fusion 
of the basisphenoid and basioccipital bones 
(Williams & Genoways 1981). The adults 
were classified into two categories based on 
the degree of fusion and ossification. The 
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skulls of older adults continue to grow, caus- 
ing shape of the skull to become more an- 
gular and massive (Russell 1968); this is 
particularly noticeable in rostral, zygomatic 
arches as well as other cranial measure- 
ments. The oldest adult age class consisted 
of animals with a total obliteration of the 
suture line separating the basisphenoid and 
basioccipital bones. The number of animals 
in this oldest age class varied among pop- 
ulations, with most being males. The largest 
sample sizes were in the younger ofthe adult 
age classes, and these were used in the mor- 
phological analysis. 

Univariate analysis of cranial characters 
(mean, range, standard error, coefficient of 
variation) was carried out using the UNI- 
VARIATE procedure of SAS (SAS Institute 
1988a, 1988b) for each sex. The relation- 
ships among the taxa were assessed using 
multiple analysis of variance (GLM pro- 
cedure), and Tukey's studentized range test 
(TUKEY's option of GLM) was used to 
identify maximally nonsignificant subsets. 
Principal component analysis (PRIN pro- 
cedure) using a correlation matrix of char- 
acters was used to identify the source of 
variation among the characters. Mean ei- 
genvectors were computed for each taxon 
and the first two components plotted. Ca- 
nonical discriminant analysis (CANDISC 
procedure) was used to compute canonical 
variates for multivariate analysis of varia- 
tion among taxa. The centroid for each tax- 
on was plotted on the first two canonical 
variates and the 9596 confidence ellipse was 
computed using the method described by 
Owen & Chmielewski (1985). Each individ- 
ual was plotted by its first two canonical 
variates in order to visualize the dispersion 
of individuals in relationship to the cen- 
troids. The contribution of each cranial 
variable used in the canonical discriminant 
analysis was determined as a percentage, us- 
ing the procedure described by Schmidly & 
Hendricks (1976). 

Phenetic relationships were assessed us- 
ing correlation and distance matrices 


generated from the character data using NT- 
SYS (Rohlf 1988, version 1.50), with clus- 
tering using UPGMA (unweighted pair 
group method using arithmetic averages). A 
minimum spanning tree also was computed 
and taxa plotted by the first three canonical 
vectors were connected. 


Results 


The diploid number of 16 specimens rep- 
resenting the new populations in Medina, 
Zavala, and Uvalde counties is 2N = 70, 
and the fundamental number is 68. The 
X-chromosome is a large acrocentric and 
the Y-chromosome is medium-sized and 
acrocentric. The morphology of the chro- 
mosomes is indistinguishable from those 
previously reported for both G. bursarius 
and G. texensis (Baker et al. 1973, Hart 
1978). 

Four of the presumptive loci were mono- 
morphic within and among all three taxa 
(Table 1). Two fixed differences were ob- 
served (PEP-1 and SOD-1) between G. bur- 
sarius and both G. texensis and the new 
taxon. No fixed differences were observed 
between С. fexensis and the new popula- 
tion. The mean heterozygosities for the new 
taxon, G. texensis, and G. bursarius are low, 
0.037, 0.041, and 0.017, respectively. The 
percentage of the loci that are polymorphic 
varied from a low in G. texensis of 27.8%, 
to 38.9% in the new taxon, and 44.4% in G. 
bursarius. Rogers’ genetic similarities were 
high when comparing С. texensis to the new 
taxon (0.915), whereas the similarities of 
these two taxa to G. bursarius were much 
lower, 0.690 and 0.648, respectively. The 
genetic distance was low when comparing 
the newly discovered populations and G. 
texensis (0.085), whereas these two taxa were 
both more distantly related to G. bursarius, 
0.352 and 0.310, respectively. 

Coefficients of variation produced in the 
univariate analysis of the nongeographic 
variation showed exceedingly high varia- 
tion in the external characters (TL, T, HF) 
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for both sexes, and were excluded from fur- 
ther morphometric analyses. The new taxon 
had the lowest measurements of the taxa 
used in this study in 22 of the 28 skull char- 
acters measured in the two sexes (Table 2). 
In the six incidences (female — BL, DIAS, 
IOC; male— СВІ, DIAS, IOC) where this 
population did not have the smallest mean 
value, one of the populations of G. texensis 
was the smallest. Analysis of geographic 
variation in individual characters using a 
single classification MANOVA test showed 
significant differences in CBL, BL, PL, 
DIAS, ZB, МВ, SB, RB, IOC, МЗ and М1 
among the populations of females; PFD and 
LN showed no significant differences in this 
analysis. Males had fewer characters dis- 
playing significant differences among the 
taxa, and these were limited to LN, MB, SB, 
RB, IOC, M3 and М1. Further analysis of 
these data using Tukey's standardized range 
test identified much of these differences to 
be attributed to relationships between only 
two or three of the taxa, especially regarding 
G. texensis, relating it to G. attwateri and 
G. personatus. Six characters displayed sig- 
nificant differences that involved all taxa in 
both males and females (Table 3): MB, RB, 
SB, M3, М1 and IOC. The new taxon was 
not significantly different from the other G. 
texensis populations in MB, SB, M1 and 
IOC in either females or males. Rostral 
breadth (RB) deviates from this clustering 
pattern in both sexes in that the new taxon's 
measurements are much smaller than those 
of G. texensis. 

The first three principal components de- 
scribe 76.7% and 85.396 of the variation 
Observed in females and males, respective- 
ly. The eigenvectors of component I were 
all positive and range from 0.110 to 0.355. 
Skull length measurements (CBL, BL, PL) 
and mastoid breadth (MB) account for most 
of the variation observed. Zygomatic 
breadth (ZB) accounts for additional vari- 
ation in males. Components II and III also 
are influenced by the variation in condy- 
lobasal length (CBL) and basal length (BL) 


. PGM-1 


Table 1.— Alleles present in three taxa. Lowercase 
letters denote alleles appearing in frequencies less than 
5%, while uppercase letters represent occurrences greater 
than 596. 
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in both sexes; however, squamosal breadth, 
mastoid breadth, and prefrontal depth are 
important characters in females. When 
plotted by the mean values for the first two 
principal components, the relationships of 
the taxa are identical to those seen for a 
similar analysis described below using the 
canonical variates, and thus are not pre- 
sented. 

There were significant differences (Ho- 
telling-Lawley's Trace: females P « F « 
0.0001; males P « F « 0.0001) in mor- 
phology among the taxa in an analysis of 
geographic variation as detected by the 
MANOVA test. The first canonical variate 
accounts for 6796 and 57% of the phenetic 
variation in females and males, respective- 
ly. The second and third variates account 
for 1796 and 996, respectively, in females; 
and 2896 and 896 in males. Three clusters 
are evident in the two dimensional plots of 
the first two canonical variates of both sexes 
(Fig. 2). Females of both samples of G. att- 
wateri cluster with the G. personatus. The 
two populations of С. texensis also cluster 
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together, with the centroid of each popu- 
lation well within the 9596 confidence ellipse 
of the other. The centroid of the new taxon 
is well outside ofthe 95% confidence ellipses 
ofall other taxa and its ellipse includes only 
the two centroids of G. texensis at the ex- 
tremes of the confidence interval (Fig. 2). It 
should be noted here that the new taxon had 
the smallest sample size (Fig. 2), and that 
this is directly observed in the larger ellipse 
resulting from the high F-value component 
of the Owen & Chmielewski (1985) for- 
mulation of the 9596 ellipse. 

Males do not show as clear a demarcation 
among the clusters as do the females al- 
though an identical pattern is evident. The 
populations of Geomys attwateri cluster to- 
gether and overlap G. personatus. The two 
populations of G. texensis form overlapping 
clusters with the centroid of the new taxon. 
The analysis of males 1s affected by smaller 
sample sizes and the increased variation as- 
sociated with the continued growth of their 
skulls throughout adult life. 

Those skull characters representing length 
(CBL, BL) and width (ZB, MB) account for 
most of the variation in canonical variate I 
of males (Table 4). Basal length (BL), palatal 
length (PL), and mastoid breadth (MB) pro- 
vide most ofthe variation seen in the second 
variate. The third variate is also highly 
weighted to length and width variables. A 
similar pattern is seen in females, but this 
length and width variation 1s seen to be lim- 
ited to condylobasal length (CBL) and mas- 
toid breadth (MB) in the first variate. Over- 
all skull length accounts for the variation 
seen in the second variate, whereas palatal 
length (PL) and squamosal breadth (SB) ac- 
count for most of the variation seen in the 
third variate. 

The phenograms constructed using the 
correlation matrices are identical between 
the sexes, and the cophenetic correlation 
values are 8796 for females and 7696 for 
males (Fig. 3). The phenograms based on 
distance matrices differ between both sexes 
and the phenograms derived from the cor- 
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relation matrices. The cophenetic correla- 
tion values are high for females (0.84) and 
males (0.87). The new taxon is quite distinct 
in both sexes, with the greatest distinction 
appearing in males. Although the clustering 
relationships vary among the analyses and 
sexes, the new taxon is always distinct when 
comparing the branch lengths to those of 
the other species. А similar relationship is 
seen in the branching pattern resulting from 
the minimum-spanning analysis (Fig. 3). 
Females of the new taxon join С. texensis, 
but G. attwateri joins G. personatus. Males 
also show a close relationship among the 
two populations of G. texensis and the new 
taxon. 


Discussion 


The newly discovered populations of 
pocket gophers at first were suspected to be 
only range extensions of taxa known to oc- 
cur in the region. However, analysis of chro- 
mosomal morphology quickly identified 
them as being quite different from any taxon 
occurring in southern Texas, and more 
closely related to gophers in the G. bursarius 
complex to the north. They share a large 
and distinctive acrocentric X-chromosome 
and 68 acrocentric autosomes with G. tex- 
ensis and race D of G. bursarius major (Ba- 
ker et al. 1973, Hart 1978). They contrast 
with Geomys personatus and G. attwateri, 
both of which have a large, subtelocentric 
X-chromosome. Furthermore, G. persona- 
tus and G. attwateri are reported to share a 
diploid number of 70 in southern Texas and 
to possess a small metacentric autosome not 
seen in the new taxon. 

In order to ascertain the relationship of 
the new taxon to С. fexensis and С. bur- 
sarius, and to determine its taxonomic af- 
filiation, we used starch gel electrophoresis 
to assay biochemical variation in 18 loci 
coding for structural proteins. Two fixed dif- 
ferences were observed between G. bursari- 
us and both G. texensis and the new taxon. 
No fixed differences, however, were ob- 
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Table 3.—Six characters that separate taxa when analyzed in single classification MANOVA. Taxa are grouped 
in nonsignificant subsets (Tukey's studentized range test) represented by the horizontal lines. Age class 4 males 
and females: characters by locality. 


n mean min max 
Males 
MB 
G. t. bakeri 5 21.96 21.4 DOR 
G. texensis 4 22.65 20.8 23.9 
G. texensis 7 23.25 21.9 23.9 
G. personatus 13 23.98 223 26.7 
G. attwateri 11 24.06 22.6 25.7 
G. attwateri 5 24.11 28.5) 24.7 
SB 
G. t. bakeri 5 16.94 16.6 17.4 
G. texensis 7 17.52 16.9 18.0 
G. texensis 4 17.63 16.8 18.5 
G. attwateri 11 18.22 17.4 19.5 
С. attwateri 5 18.63 18.1 19.5 
G. personatus 13 18.69 17.5 19.8 
RB 
G. t. bakeri 5 9.19 8.9 9.5 
G. personatus 13 0217 8.9 10.7 
С. attwateri 11 9.43 9.1 9.9 
G. attwateri 5 9.74 8.9 10.4 
G. texensis 4 9.95 9.8 10.2 | 
G. texensis 7 10.17 9.3 10.9 
M3 
G. t. bakeri 5 6.90 6.8 23] | 
G. texensis 4 7.48 p 7.8 
G. texensis 7 7.61 7.3 8.0 
G. attwateri 11 7.74 7.1 8.1 
G. personatus 13 7.80 7.3 8.8 
G. attwateri 5 7.92 7.6 8.5 
MI 
G. t. bakeri 5 8.49 8.2 8.8 
G. texensis 4 8.63 8.1 9.0 
G. texensis 7 8.91 8.6 9.6 
G. attwateri 11 9.36 9.1 9.8 
G. attwateri 5 9.36 9.1 9.6 
G. personatus 13 9.43 8.4 10.6 
IOC 
G. texensis 4 5.87 5.5 6.3 
G. texensis 7 5.88 5.6 6.2 
С. t. bakeri 5 5.96 5.7 6.4 
G. personatus 13 6.18 5.8) 6.6 
G. attwateri 11 6.40 6.0 6.7 
С. attwateri 5 6.53 6.1 6.9 


served between С. texensis and the new tax- geomyids from central Texas. Genetic sim- 
on. Block & Zimmerman (1991) identified  ilarities determined by them for С. bursari- 
fixed differences between С. bursarius and us and populations of С. texensis ranged 
G. texensis in a study involving species of from 0.607 to 0.648, whereas genetic sim- 
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Table 3.— Continued. 
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n mean min max 
Females 
MB 
G. t. bakeri 6 21.14 20.5 21.9 
G. texensis 15 21.36 19.9 22.8 
G. texensis 9 Meth 20.8 23.0 
G. personatus 13 22В98 207 24.0 
С. attwateri 18 23.30 21.6 24.8 
G. attwateri 3 23.30 22.8 24.1 
SB 
G. t. bakeri 6 16.65 16.3 17.1 
G. texensis 9 16.85 16.5 17.3 
G. texensis 15 16.92 16.0 17.9 
G. attwateri 18 17.97 16.9 18.7 
G. attwateri 3 18.01 17.7 18.3 
С. personatus 13 18.21 17.3 18.9 
RB 
G. t. bakeri 6 8.40 8.0 8.8 
G. attwateri 18 8.56 8.2 9.7 
G. personatus 13 8.87 8.1 9,5 
G. texensis 15 9.15 8.5 9.6 
G. texensis 9 9.37 9.0 10.1 
G. attwateri 3 9.60 9.0 10.1 
M3 
G. t. bakeri 6 6.88 6.5 7al 
G. texensis 15 7.36 7.0 7.8 
G. attwateri 3 7.50 7.4 7.6 
G. attwateri 18 7-51 7.0 8.0 
G. texensis 9 7.52 7.1 8.1 
G. personatus 13 7.72 p» 8.2 
MI 
G. t. bakeri 6 8.25 7.8 8.5 
G. texensis 15 8.51 7.9 9.3 
G. texensis 9 8.66 8.1 9.0 
G. attwateri 18 8.99 8.5 9.4 
G. personatus 13 9.08 8.5 9.5 
G. attwateri 3 9.21 9.0 DAS 
IOC 
G. texensis 9 5.76 DES 5.9 
G. t. bakeri 6 5.83 5.5 (285 
G. texensis 15 5.86 5.4 6.3 
G. personatus 13 6.07 S.8 6.6 
G. attwateri 18 6.28 5.6 6.7 
G. attwateri 3 6.49 6.2 6.9 


ilarities produced from intraspecific com- 
parisons involving G. texensis were high 
(from 0.931 to 0.937). The new taxon has 
a high genetic similarity when compared to 
G. texensis (0.915). Similar comparison to 


G. bursarius yields a much lower value 
(0.648). Similarity values of 0.9 are com- 
parable to previous studies of geomyids 
when making intraspecific comparisons of 
populations, and the lower genetic similar- 
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ity values, seen when comparing the new 
taxon to G. bursarius, are within the range 
associated with interspecific comparisons 
(Block & Zimmerman 1991, Dowler 1982). 

Relationship between G. texensis and the 
new taxon is supported further by the mor- 
phometric analysis. In an analysis of the 
cranial characters that appear to separate G. 
personatus and G. attwateri from G. tex- 
ensis, the new taxon was associated with G. 
texensis in mastoid breadth, squamosal 
breadth, interorbital constriction, and length 


of molar toothrow. A relationship to G. tex- 
ensis also is seen in the two phenetic clus- 
tering results. The phenograms based on the 
correlation matrices describe the new taxon 
as similar to the two G. texensis popula- 
tions. Minimum-spanning analysis yields 
similar relationships with a population of 
G. texensis serving as the branching neigh- 
bor to the new taxon. 

These data support a proposal that the 
new taxon is closely related to G. texensis. 
The question then arises, do these new pop- 
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Table 4.— Coefficients for canonical variates and the percent influence of each variable for the variates. 


Canonical Canonical Canonical 
Variate I Variate II Variate III 
Variable Variable Variable 
Character coeff. % coeff. % coeff. % 
Males 
CRL — 1.8433 28.51 —0.3363 3.36 —4.3026 41.24 
BL 0.6980 10.14 —4.5368 42.57 1.8989 17.09 
PL —0.6005 6.09 2.1302 13.95 207230 17.11 
PFD —0.9800 5.87 —0.6656 2.58 0.9685 3.60 
DIAS 1.1349 6.15 2.1888 7.67 —0.5400 1.81 
ZB —1.1519 11.26 0.6130 3.87 0.5727 3.47 
МВ 2.0305 17.87 2.2454 12.77 —2.0193 11.01 
SB 0.8427 5.67 — 1.8478 8.03 0.3642 1.52 
RB —0.5925 2.10 0.0785 0.18 —0.1217 0.27 
IOC —0.6612 1.52 0.8264 1.23 0.6400 0.91 
M3 0.8382 2.38 — 1.2829 2.35 0.6570 1.16 
MI 0.7154 2.44 0.6592 1.45 0.3881 0.82 
Females 

CBL — 1.9191 35.85 0.7594 15.53 — 0.0083 0.28 
BL 0.4529 7.94 — 1.7802 34.17 —0.2815 8.76 
PL 0.5315 6.44 0.9833 13.05 1.3095 28.18 
PFD —0.4949 3.63 0.1148 0.92 —0.3503 4.56 
DIAS 0.6719 4.25 — 1.1497 7.96 —0.7237 8.13 
ZB 0.1029 1.18 0.3581 4.49 0.5166 10.49 
MB 2.4876 26.45 —0.5207 6.06 0.5198 9.81 
SB —0.3687 3.06 0.9749 8.87 — 1.3019 19.20 
RB — 1.0596 4.56 0.7393 3.48 0.6176 4.71 
IOC 0.7276 2.09 — 0.2932 0.92 0.7603 3.87 
M3 —0.4172 1.50 0.9869 3.88 — 0.0849 0.54 
М1 0.7203 3.05 0.1459 0.68 —0.1930 1.45 


ulations, 1solated and 120 km distant from 
the previously known distribution of G. tex- 
ensis, constitute a new subspecies? An anal- 
ysis of cranial morphology identifies the new 
populations as being quite distinct from G. 
texensis. They are smaller in 22 of 28 com- 
parisons involving both sexes. Canonical 
discriminant analysis identifies measure- 
ments reflecting the length of the skull (BL, 
PL) and mastoid breadth accounting for 
much of the variation separating the taxa. 
The centroids produced from plotting these 
taxa along the first two canonical variates 
are widely separated, although the ellipses 
are overlapping. Separation 15 more clearly 
seen in analysis of females, where the two 
populations of С. texensis are closely as- 
sociated with each other, and each is well 


within the ellipse of the other. The new tax- 
on has a centroid separated from these two 
primarily along the axis of the second ca- 
nonical variate. The centroid 15 well outside 
the range of the ellipses of the G. texensis 
populations although the two centroids of 
G. texensis occur within the 9596 ellipses of 
the new taxon. This distinction is less clear 
in males where the ellipses of both taxa in- 
clude the centroids of the other. The new 
taxon 15 also seen to be quite distinct from 
the two populations of G. texensis in the 
phenogram produced from the clustering of 
the correlation matrices. 

The populations of ©. texensis represent 
two previously recognized subspecies, G. b. 
texensis and G. b. llanensis. Honeycutt & 
Schmidly (1979) identified primarily size- 


18 PROCEEDINGS OF THE BIOLOGICAL SOCIETY OF WASHINGTON 


o 
о 
CAN il 


Females 
—2.0 
2.0 
s td 
= 
0.0 
5 


6. attwateri 


—2.0 0.0 2.0 4.0 


Females 
0.99 
‚к —+ 
Correlation Matrix 

. p. streckeri 


. attwateri 


. texensis 


G 
G 
G. attwateri 
G 
G. texensis 
G 


. t. bakeri 


Females 

DN ner 

Distance G. p. streckeri 
С. texensis 
G. texensis 
G. t. bakeri 
G. attwateri 
б. attwateri 

Fig. 3. 


CAN Ill 


G. p. etreckeri 


Males 


0.99 
—— 


Correlation Matrix 


G. p. streckeri 
G. attwateri 
G. attwateri 
G. texensis 
G. texensis 
G. t. bakeri 
Males 
12778) 1.00 0.25 
Distance 


G. p. streckeri 
. texensis 
. texensis 


. attwateri 


oO o о 9 


. attwateri 


G. t. bakeri 


Plot of the centroids of taxa along their first three canonical variates. The centroids are connected 


by the branching order of the minmum-spanning analysis produced from the variance-covariance matrices of 
measurements of cranial characters. UPGMA phenograms produced from correlation and distance matrices of 


the cranial measurements. 


related differences between these taxa. When 
compared to these populations of С. tex- 
ensis, the new populations seem to have dif- 
ferences in cranial features that have re- 
sulted in alteration of both skull size and 
shape. Cranial morphology has been pro- 


posed to be greatly influenced by both soil 
composition and texture (Hendricksen 
1972, Smith & Patton 1988, Wilkins & 
Swearingen 1990). Wilkins & Swearingen 
(1990) noted an increase in the mean values 
of all cranial characters in populations of G. 
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personatus in sandy soils when compared 
to other soil types. This difference also ex- 
tended into a multivariate analysis which 
effectively separated populations inhabiting 
fine sandy loams, loamy fine sand, and fine 
sand, with an increase in size from smaller 
to larger along that axis. Geomys texensis 
occurs in sandy-loam soils in the central 
basin region of the Edwards Plateau. These 
are porous, well drained soils. The new tax- 
on from South Texas inhabits a loam sand, 
Atco soil that is denser and less friable. Go- 
phers in the heavier and denser soil are 
smaller and have cranial changes which al- 
ter the skull shape, and conform to the mod- 
els described by Smith & Patton (1988) and 
Wilkins & Swearingen (1990). 

Geomys texensis was previously known 
only from the central basin of the Edwards 
Plateau, and isolated there by barriers of 
stony to gravelly clay, clay, and shallow 
loamy soils (Godfrey et al. 1973). Block & 
Zimmerman (1991) described a scenario in- 
volving a late Wisconsin to Holocene en- 
vironmental change that isolated G. tex- 
ensis as the warmer and drier conditions 
approximately 9000 B.P. accelerated ero- 
sion. The presence of G. texensis along the 
southern edge of the Edwards Plateau is 
plausible within the framework of this hy- 
pothesis. The distribution of G. texensis once 
could have been more widespread in south- 
central Texas, and probably extended fur- 
ther into southern Texas. A cooler climate 
and deeper soils would have allowed G. tex- 
ensis in the central basin to contact popu- 
lations south of the Edwards Plateau across 
the upper reaches of the Llano River drain- 
ages. Subsequent xeric conditions merely 
separated these populations, isolating those 
to the south and ultimately restricting them 
to the smaller pockets of suitable soils south 
of the newly created indurate soils of the 
plateau. Fossil remains of geomyids from 
cave deposits on the Edwards Plateau 
(Dalquest & Kilpatrick 1973) provide evi- 
dence for a wider distribution, with geo- 
myids ranging over at least the eastern por- 
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tions of the plateau from 10,000 to 4000 
B.P. 

The new subspecies described below 
demonstrates close affinities to G. texensis 
in cranial and chromosome morphologies, 
and in biochemical variation. Analyses of 
cranial morphology indicate that these pop- 
ulations resemble G. texensis, but appear to 
differ significantly in having size- and shape- 
related changes. These cranial differences 
alone could indicate a species distinction 
but the presence of a low level of genetic 
differentiation leads us to be more conser- 
vative. We conclude that, based on the ob- 
served variation and the extreme spatial 
separation, this new taxon 1s related to G. 
texensis and is an isolated, relictual popu- 
lation of this species, forming a distinctive 
subspecies. 


Geomys texensis bakeri, new subspecies 


Holotype. — Adult male, skin, skull, and 
body skeleton, no. 52310, Texas Coopera- 
tive Wildlife Collections (TCWC); from 1 
mi E D'Hanis, Medina Co., Texas; obtained 
on 3 Jan 1987 by R. M. Pitts, original no. 
1998. 

Distribution. — Two isolated populations 
have been found along separated drainages 
of the Frio River. One population occurs 
along the Sabinal and Frio rivers in Uvalde 
and Zavala counties. А second population 
in Medina County is restricted to soils along 
Seco and Parker creeks, tributaries of the 
Frio River. Both populations are associated 
with nearly level Atco soil (Stevens & Rich- 
mond 1976, Dittmar et al. 1977), which has 
a patchy distribution in this region. The soil 
is well drained and consists of sandy surface 
layers with loam extending to as deep as 2 
m. This soil is associated with stream ter- 
races formed by the drainage systems in each 
locality. These two populations are widely 
separated (40 km); however, there may be 
additional populations along Seco and Par- 
ker creeks as they flow southward to join 
the Frio River. The nearest geomyids are 
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G. attwateri (Medina County), G. persona- 
tus streckeri (Dimmit County), and G. t. tex- 
ensis (Kimble County). 

Description. —SSmall size which is es- 
pecially evident in the measurements of 
body length, skull length (CBL, BL), nasal 
length (LN), and skull breadth (ZB, MB, 
SB). Pelage coloration is russet brown on 
the dorsum, and grades to a paler color along 
the sides. The basal portions of the hairs are 
gray. A dark dorsal stripe extends from the 
head to the rump. The ventral surface is 
white with gray coloration on basal parts of 
the hairs. The tail is sparsely haired and 
consists of a mixture of brown and white 
hairs. The feet are white haired. Subadult 
pelage is a tawny brown, whereas adult pel- 
age is a darker, richer brown, and appears 
glossy and more reflective. 

Pelage color appears paler in gophers col- 
lected in the more sandy surface soil along 
the Frio River in Uvalde and Zavala coun- 
ties then in those from Medina County. This 
difference in pelage coloration correlates 
with the much paler color of the substrate 
at this locality than the darker color of the 
loamy soil at the Medina locality. 

Karyotypic features. —'The diploid num- 
ber is 70 and the fundamental number is 
68. The X-chromosome is a large acrocen- 
tric and the Y-chromosome is a medium- 
sized acrocentric. 

Measurements. — Measurements (in mil- 
limeters), as described in Williams & Geno- 
ways (1977), from 12 adult individuals from 
the two populations are listed in Table 4. 
Measurements of the holotype (TCWC 
52310) are: total length, 227; length of tail, 
66; length of hind foot, 27; length of ear, 6; 
greatest length of skull, 40.8; condylobasal 
length, 40.4; basal length, 38.3; palatal 
length, 26.3; prefrontal depth, 15.5; length 
of nasals, 13.3; length of diastema, 14.3; 
zygomatic breadth, 25.6; mastoid breadth, 
22.4; squamosal breadth, 16.9; rostral 
breadth, 9.4; interorbital constriction, 6.3; 
breadth across third molars, 7.1; length of 
maxillary toothrow, 8.5. 


Comparisons. — Cranial measurements of 
individuals of G. t. bakeri are smaller in size 
than those of G. attwateri, G. p. streckeri 
and G. bursarius major (Baker & Genoways 
1975). This is especially evident in the mea- 
surements reflecting the length (CBL, BL) 
and breadth (ZB, MB, SB) of the skull. Geo- 
mys b. major is a larger gopher in all com- 
parative external and cranial measurements 
as seen when comparing the measurements 
herein to those given by Baker & Genoways 
(1975). Geomys attwateri 15 the taxon near- 
est geographically to G. t. bakeri. It has ex- 
ternal measurements comparable to those 
of G. t. bakeri, but the skull is longer (CBL, 
BL) and wider. Significant differences are 
detectable in mastoid breadth, squamosal 
breadth, and interorbital constriction (С. t. 
bakeri males MB = 21.9, SB = 16.9, IOC 
= 5.9; females MB = 21.1, SB = 16.6, ТОС 
= 5.8; С. attwateri males MB = 24.0, SB = 
18.6, ТОС = 6.4; females MB = 23.3, SB = 
18.0, ТОС = 6.3). Geomys attwateri has a 
paler pelage that 15 a buffy tan in color, and 
the pelage has a uniform, nonglossy ap- 
pearance. Other populations of G. attwateri 
in the eastern part of southern Texas have 
pelage colors that are similar to that of G. 
t. bakeri. 

The karyotypes of Geomys attwateri and 
G. personatus differ from that of G. t. bakeri 
in that they have a large, submetacentric 
X-chromosome and a small, metacentric 
autosome (Davis et al. 1971, Honeycutt & 
Schmidly 1979, Tucker & Schmidly 1981). 
Geomys texensis (Honeycutt & Schmidly 
1979) and races of Geomys bursarius (Baker 
et al. 1973) have an identical diploid and 
fundamental number, and appear identical 
when standard karyotypes are compared. 

Etymology. — The subspecific name is a 
patronym and is selected to honor Dr. Rob- 
ert J. Baker for his many contributions to 
mammalogy, particularly to the systematics 
and evolution of Geomys, as well as his 
overall research program, activity in pro- 
fessional societies, and involvement in 
graduate education. 
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Appendix 


Specimens examined. — Two hundred ninety seven 
specimens were used from the following collections: 
American Museum of Natural History (AMNH), Mid- 
western University (MSU), Tarleton State University 
(TSU), Texas Cooperative Wildlife Collections 
(TCWC), Texas Natural History Collection (TNHC), 
Texas Tech University (TTU), National Museum of 
Natural History (USNM). 

Geomys attwateri. — (85). Texas: Atascosa Co.: 6 mi 
W Campbellton, 1(TNHC); 2 mi NW Campbellton, 
3(TCWC); Е mi E Lytle, 16(TCWC); 7 mi E Lytle, 
16(TNHC); 2.4 mi SE Lytle, 3(TCWO); 7 mi SE Lytle, 
8(TCWC). Bexar Co.: 15 mi SE San Antonio, 1 (TSU). 
Frio Co.: 1 mi N Moore, 3(TCWO); Pearsall city limits, 
2(TCWO); 2.25 mi S, 1 mi E Pearsall, 1(TCWC); Mc- 
Coy, 1(TNHC); 0.3 ті E McCoy, I(TNHC); 2 mi N 
Pleasanton, 7(ТМНО. Medina Co.: 5 mi W Devine, 
4(TCWC); 7.2 mi E Yancy, l(TCWC), 1(TSU). Wilson 
Co.: 11 mi NW Floresville on HWY 181, 6(TCWC); 
4 mi W Floresville, 1(TNHC); 1 mi W Floresville, 
2(TCWO); 34 mi S, 2⁄2 mi E Floresville, 4(TCWC); 3.2 
mi NW Poth, 1(TNHO); 3.6 mi SSE Poth, 1(ТМНС); 
5.4 mi W San Antonio River, between Floresville and 
Pleasanton, 1(TNHC). 

Geomys personatus streckeri. — (68). Texas: Dimmit 
Co.: 13 mi М Carrizo Springs on HWY 277, 3(1TU); 
13 mi NE Carrizo Springs, 3(TTU); 13 mi NE Carrizo 
Springs on US HWY 277, 5(TTU); 4 mi N, 9 mi W 
Carrizo Springs, 1(TCWC); 1 mi S Carrizo Springs, 
700 ft, 2(TCWC); 15 mi S, 11 mi W Carrizo Springs, 
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5(TCWC); 1!^ mi E Carrizo Springs, 6(TCWC); E Car- 
rizo Springs, 3(TTU); Carrizo Springs, 30(TNHC), 
1 1(TCWO); near Carrizo Springs on HWY 277, 5(TTU); 
1.0 mi SW Carrizo Springs, HWY 186, 4(TCWC); 2 
mi S Carrizo Springs, low water crossing Dentonio Rd, 
I(TCWC). 

Geomys texensis bakeri. —(35). Texas: Medina Co.: 
1 mi E D'Hanis, 5(TCWC); D'Hanis, 5(TCWC); 5% 
mi W Hondo, 6(TCWC); 6.2 mi W Hondo, 4(TCWC). 
Uvalde Co.: 16 mi S Sabinal on FM 187, 6(TCWC); 
17 mi S Sabinal on FM 187, 4(TCWC). Zavala Co.: 
18 mi $ Sabinal on FM 187, 1(TCWOC); Ио mi S Uvalde 
County line, 4(TCWC). 

Geomys texensis texensis. —(108). Texas: Gillespie 
Co.: 13 mi N Fredericksburg, 1(TSU); 1 mi N Fred- 
ericksburg, 2(TCWC); 0.5 mi М Fredericksburg, 
I(TCWC); 9 mi W Fredericksburg, 1(TNHC). Kimble 
Co.: Junction, 3(TCWC). Llano Co.: 2.6 mi N, 1.8 mi 
E Castell, S(TTU); 6.4 mi E Castell, 1(TCWC); Castell, 
I(TTU); 6.4 mi E Castell 1(TCWC); 1 mi E Castell 
I(TCWC); 1.2 mi W Castell on FM 152, 1(TCWC); 8 
mi S, 0.9 mi W Kingsland, 4(TTU); 9.2 mi S, 1.1 mi 
E Kingsland, 1(TTU); 10 mi S, 1.8 mi E Kingsland, 
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2(TTU); 2.9 mi NW Llano on HWY 71, 2(TTU); 0.2 
mi М, 8.7 mi W Llano, 3(TTU); Drier Cr at Lone 
Grove, 7 mi W Llano, IlO(TCWC); 4 mi W Llano 
I(TCWC), Llano, 1(TCWC); 0.2 mi E Llano, 1(TCWC); 
1 mi E Llano, 2(TCWC); 2 mi E Llano, 4(TCWOQ); 7.2 
mi E Llano, 1(TCWC); Oatman Cr, 3 mi S Llano, 
6(TCWC); 51.6 mi W Austin on HWY 71, 2(TTU); 3 
mi S Jct FM 268 and HWY 29 on 29, 2(TTU); 9.3 mi 
N Jct Texas 29, Texas 16 on Texas 16, 1(TTU). Mason 
Co.: 3.1 mi E Art, 1(TCWO); Art city limits, 3(TCWC); 
12 mi М Mason, 6(MWU); 3.6 mi М, 1.5 mi W Mason, 
I(TTU); 1 mi М, 1.1 mi W Mason, 4(TTU); 12 mi W 
Mason, 2M WU); 9.4 mi W Mason on US 377, 3(TTU); 
Mason, I(TCWC); 1 mi E Mason, 6(TCWC); 6.5 mi 
E Mason on Texas 29, I(TTU); 2.0 m S, 2.7 m W 
Castell, 1(TTU); 5 mi S Mason, J(MWU). in Mason 
Co.: 3% mi W Castell, 3(TCWC); 0.3 mi S, 1.5 mi W 
Castell, 3(TTU); 0.3 mi S, 0.9 mi W Castell, 1(TTU); 
0.3 mi S, 0.8 mi W Castell, 1(TTU); 2.6 mi S, 3.0 mi 
W Castell, I(TTU); 0.7 mi S, 2.1 mi W Castell, 2(TTU); 
1.0 mi S, 2.3 mi W Castell, 2(TTU); 11 mi NE London, 
HWY 377, 2(TCWC); 13 mi NE London, HWY 377, 
2(TCWC). 


